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ABSTRACT 

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive 
constraints on its dark matter distribution. Though the data are unable to distinguish between models with small 
cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. 
Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its 
core radius to be no larger than ^ 700 pc. We derive an upper limit of rcore ^ 300 pc by demanding that the 
central phase-space density of Fornax not exceed that directly inferred from the rotation curves of low-mass 
spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the 
core to be quite small in order to avoid conservative limits from the Lya forest power spectrum, rcoie 85 pc. 
We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can 
be explained by the presence of a large ^ L5 kpc core. A self-consistent core of this size would be drastically 
inconsistent with the expectations of standard warm or cold dark matter models, and would also require an 
unreasonably massive dark matter halo, with Vm-dx — 200 kms"' . 

Subject headings: cosmology: dark matter, cosmology: observations, cosmology: theory, galaxies: kinematics 
and dynamics, galaxies: structure, galaxies: formation, dwarfs: galaxies 



1. INTRODUCTION 

The central density distribution in many dark matter dom- 
inated galaxies appears to be lower than what is expected 
in standard A + Cold Dark Matter (ACDM) models (Moore 
1994; Flores & Primack 1994; de Blok et al. 2001; Zentner & 
Bullock 2002; Kuhlen et al. 2005; Simon et al. 2005; Zackiis- 
son et al. 2006). This discrepancy, along with other potential 
difficulties (e.g. Klypin et al. 1999; Moore et al. 1999), can 
be ameliorated by considering alternative dark matter models 
(e.g. Zentner & Bullock 2003, and references therein). One 
particularly intriguing scenario is Warm Dark Matter (WDM), 
which differs from CDM in that the dark matter has a non- 
negligible velocity dispersion (Bond et al. 1980; Blumenthal 
et al. 1982; Pagels & Primack 1982; Dodelson & Widrow 
1994; Hogan & Dalcanton 2000; Abazajian 2006a). Gener- 
ally, these WDM models come with two distinct observational 
signatures. First, there is a reduction in clustering and a delay 
in collapse times for structures on linear scales approximately 
smaller than the free-streaming scale of the WDM particle. 
This effect is currently probed by measurements of the Lya 
forest power spectrum and this places strong constraints on 
the mass of the WDM particle (Viel et al. 2005; Abazajian 
2006b; Seljak et al. 2006). Second, the velocity dispersion 
of WDM imposes an upper limit on the phase-space density, 
defined as 2 = p/'^^^ where p is the density and <j is the ve- 
locity dispersion of the dark matter The upper limit on Q 
implies that the density profiles of WDM halos must saturate 
to form a constant density core at an inner radius defined by 
Q < 6max (Tremaine & Gunn 1979). This is in stark con- 
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trast to CDM predictions that give dark halo density profiles 
that rise steeply towards the center, p oc r'', with 7 ~ 1, (e.g. 
Dubinski & Carlberg 1991; Navarro et al. 1996; Moore et al. 
1999; Klypin et al. 2001; Navarro et al. 2004; Diemand et al. 
2005) and have divergent phase-space profiles Q cx r"' -' (Tay- 
lor & Navarro 2001). 

Dwarf spheroidal galaxies (dSph's) are diffuse, low- 
luminosity systems, with a total mass believed to be domi- 
nated by their host dark matter halos (Mateo 1998). Assum- 
ing that these galaxies are in dynamical equilibrium, the stars 
act as tracers of the gravitational potential, and can be used 
as a probe of the dark matter density profile. Indeed, dSph's 
provide a unique testing ground for the nature of dark mat- 
ter, because in these low-mass systems the phase-space cores 
are the most observationally accessible (Sellwood 2000; Dal- 
canton & Hogan 2001; Alam, Bullock, & Weinberg 2002). 
Unfortunately, attempts to constrain uniquely the dark mat- 
ter profiles of dwarfs using stellar velocity dispersion data are 
plagued by degeneracies. We emphasize here that while the 
data are currently unable to settle the question of a central 
density core versus cusp, they do provide powerful constraints 
on the combination of core size (rcm-e) and halo maximum cir- 
cular velocity (Vmax)- Each solution maps directly to a con- 
straint on Q. 

The Fornax dSph is approximately 138 kpc from the Milky 
Way (e.g., Buonanno et al. 1999), in the vicinity of its orbital 
pericenter (Dinescu et al. 2004). Its stellar population is dom- 
inated by intermediate age stars (^ 4 Gyr), although Fornax 
has also both very young stars (^ 0.2-2 Gyr) and very old 
stars (> 10 Gyr, e.g., Pont et al. 2004, and references therein). 
Among the most puzzling observed properties of Fornax is 
the fact that it hosts 5 globular clusters distributed over a wide 
range of distances from the center of this galaxy. Standard ar- 
guments suggest that (in the absence of external heating) these 
systems should have sunk to the center via dynamical friction 
in much less than a Hubble time (Oh et al. 2000). Recently, 
Goerdt et al. (2006, G06 hereafter) and Sanchez-Salcedo et al. 
(2006) suggested that the presence of these clusters provides 
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Fig. 1 . — Measurements of the Fornax velocity dispersion profile compared to models for the dark matter halo. Left panel: Predictions for the case of a fixed 
''core = 1 kpc; the solid curve shows a large Vmax = 140 km/s, and the dashed curve shows a small y„,ax = 30 km/s. Right panel: The best-fitting cases for two 
distinct density profiles: the solid curve shows the fit for the a = 1.5 profile described in the text with Vmax = 28 km/s, and the long-dashed curve shows the 
best-fitting cuspy NFW profile with Vmax = 50 km/s. 



good evidence for the existence of a large, constant density 
dark matter core in the center of Fornax. Specifically, G06 ar- 
gued that the globular clusters will stall in their orbital decay 
if they encounter a constant density core. Thus the globular 
clusters should be located in a narrow range of radii from the 
center of Fornax corresponding to its halo core size. The For- 
nax globular clusters are observed to orbit at a range of pro- 
jected distances from the center of galaxy of between 0.24 and 
1 .60 kpc (Mackey & Gilmore 2003). This suggests that a core 
size rcoie ^1-5 kpc would be required to explain the observed 
distribution. Another idea, which relies on a somewhat fine- 
tuned timing argument, is that only the central-most globular 
cluster is stalled at the core, rcore ^ 240 pc, (G06; J. Read, 
private communication) and the other four globular clusters 
were simply formed at a large enough radius that they have 
not sunk to the center. 

In this paper we consider the general plausibility of a large 
core in Fornax. In section 2, we focus on the allowed dark 
halo solutions which reproduce the observed stellar kinemat- 
ics, and we present a toy model to understand the nature of the 
these solutions. In section 3, we combine these solutions with 
limits on the WDM mass from the Lya forest power spec- 
trum to place constraints on the maximal core size. We also 
discuss additional, more independent constraints on the core 
size. We conclude by discussing these results in the context 
of the observed globular cluster distribution in Fornax. 

2. DYNAMICAL CONSTRAINTS ON THE DARK HALO OF FORNAX 

For a stellar system in equilibrium and embedded in a 
spherically symmetric dark matter potential, the radial stellar 
velocity dispersion (Jrir) is given by the Jeans equation: 



r 

dr 



= -p,{r)V^(r)-2(3(r)p*'Jr- 



(1) 



Here the circular velocity is governed by the mass distribu- 
tion V^(r) = GM{r)/r and the anisotropy parameter is = 
1 - (Tg/ a^, where (Tr and ae are the radial and angular veloc- 
ity dispersions, respectively. We consider (iir) to be constant 
in radius, and explore three choices of /3 = -0.5, 0, and 0.5, 



using /3 = as our fiducial case. The appropriate value for (3 
is unknown, however this range brackets reasonable choices. 
Dissipationless CDM simulations show that (3 for dark mat- 
ter increases from ^ to 0.6 as the radius increases towards 
the virial radius of the halo (Cole & Lacey 1996; Colin et al. 
2000; Diemand et al. 2004; Faltenbacher et al. 2005; Wojtak 
et al. 2005). However, models of dSph galaxies (e.g. Mayer & 
Wadsley 2004) suggest that a more natural value for the stars 
may be negative, (3 ^ -0.5. Our main results are not sensitive 
to the choice of /3 > 0, and we note that /? < provides an 
even stronger upper limit on the core size than the /3 = case. 

We assume that the stellar component and the dark matter 
component are uncoupled, each with an independent density 
distribution, and that M{r) is dominated by the dark matter. 
We model the stellar distribution, /?*(r), as a spherically sym- 
metric King profile (King 1962) and adopt = 0.39 kpc and 
rt = 2.7 kpc for its King core and tidal radii (Mateo 1998). For 
the dark matter, we consider a density profile of the form 



Pir) = 



Po 



[l + (r/ror 



|3/q ■ 



(2) 



where po is a central core density and ro is a characteristic ra- 
dius. Note that we demand that the profile falls off at large 
radius as r"^, as seen in isolated galaxies (Prada et al. 2003). 
The parameter a controls the sharpness of the transition from 
flat to falling density with radius. In our fiducial case we 
adopt a = 1.5. This choice yields a profile that matches 
very closely that advocated by Burkert (1995) to match rota- 
tion curve data for low-mass galaxies. However, our conclu- 
sions do not change considerably if we adopt other reasonable 
choices a = 1.0 or 2.0 (see below). Note that once a is fixed, 
any two independent parameters in the mass distribution (e.g. 
ro and the maximum of the circular velocity, Vmax) define the 
density profile completely. The circular velocity curve im- 
plied by equation 2 peaks at a radius rmax = (4.4, 3.2, 2.9) ro 
for a = (1.0, 1.5,2.0). 

The log-slope of this profile, 7, gradually approaches zero 
towards the center, but a constant density "core radius" is 
never achieved at finite r. A natural way to define a core radius 
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Fig. 2. — Constraints on the core radius of Fornax as a function of central phase-space density (left) and maximum circular velocity (right) derived from the 
velocity dispersion profile. We define Tcmc as the radius where the log-slope of the density profile is 7 = —0.1 in equation 2 with a = 1.5, rcoic ~ O.lro. The 
long-dashed, solid, and short-dashed lines use /3 = 0.5, 0.0, and —0.5 respectively. Values of rcorc with Q ^ 3 X 10"^ Mq pC^ (km/s)"^ are ruled out by Lya 
forest constraints in the case of WDM with g = 2 (rcore ^ 85 pc). In a more general class of dark matter models, directly observed phase-space limits from 
low-mass spiral rotation curves, Q > 10^ Mq pc"^ (km/s)"^ demand ''core ^ 300 pc. Most generally and irrespective of the cause of the core, a Fornax halo 
with Vmax ^ 100 km s"' is disfavored from dynamical considerations. This implies i-coic ^ 700 pc. 



in this case is to use the radius where 7 reaches some small 
value, such that rcore = r^, where is defined as follows: 

r^ = ro [-3/7-1]"'/". (3) 

More general forms for p(r) may be considered, but this log- 
slope definition of core radius is always possible. In this paper 
we define the core radius to be where 7 = -0.1. This is the 
same definition used by G06. For a = 1.5 we have rcore — 
0.1 ro. 

In order to compare to measurements of the radial stellar 
velocity dispersion profile, we project equation 1 along the 
line-of-sight to the observer Performing this projection gives 

where I{R) is the surface density as determined from p^„ and 
R is the radial distance from the center of the galaxy. The line- 
of-sight stellar velocity dispersion was recently measured by 
Walker et al. (2005). The authors utilize a primary sample 
ofN = 209 stars and perform cuts to remove interloping stars, 
and present data for three removal schemes that yield = 179, 

= 185, and = 189 stars respectively. For our analysis, 
we use their = 185 sample, as shown by the points with 
with error bars in Figure 1 . Our conclusions are not strongly 
affected by this particular choice. 

The lines in Figure 1 show the predicted a-j, profiles for 
four different input dark matter density distributions com- 
pared to the data. In each case we have assumed /3 = 0. 
Each curve in the left panel assumes a = 1.5 in equation 2 
and fixes the core size to rcore = 1 kpc. The solid line shows 
the best-fit case to the data with these restrictions and has 
Vmax =140 kms"'. The dashed line, on the other hand, cor- 
responds to what might be considered a more reasonable halo 
mass, with V^ax = 30 kms"'. Clearly, in order to reproduce a 
reasonable stellar velocity dispersion profile with such a large 
core, the dark matter halo of Fornax must be quite massive. 
In order to match the data with V^nx = 30 kms"', a smaller 



core must be used. The solid line in the right panel uses 
the same profile shape with V^ax = 30 kms"', but now with 
/"core = 0.15 kpc. There is a clear degeneracy between Vmax and 
'"core - a large flat density distribution requires a deeper poten- 
tial well in order to reproduce the data. The dashed line in the 
right hand panel shows a best-fit case for a (cuspy) NFW dark 
matter profile. Note that these data do not exclude an inner 
cusp. 

Figure 2 illustrates a more general exploration of the al- 
lowed parameter space for a= 1.5. We demand that the al- 
lowed fit region obey < Xmin + ^X'^' where Ax^ = 4.6 de- 
fines the 90% confidence level for a distribution with two 
degrees of freedom. The parameters rcore and Vmax define the 
density profile and must lie within the banana-shaped con- 
tour in the right panel in order to reproduce the observed &p. 
At large Vmax the contours follow rcore Vmax and we discuss 
how this can be understood analytically below. The left panel 
shows the same best-fit parameter space, now expressed in 
terms of rcore and the central (r = 0) phase-space density of 
the dark matter, Q, implied by each solution. The different 
line types in Figure 2 correspond to different assumed val- 
ues for the velocity anisotropy: f3 = -0.5 (short-dash), (3 = 
(soHd), and /3 = 0.5 (long-dash). Note that for /3 = 0.5 solu- 
tions with rcore ^ 250 pc and Vmax 30 kms"' are disfavored. 
For P = -0.5 solutions with Vmax ^ 30 kms"' are disfavored. 

Similar results for the velocity anisotropy for Fornax were 
uncovered in previous studies (Lokas 2002; Kazantzidis et al. 
2004). In particular Lokas (2002) showed that the degeneracy 
with P and the density profile can be broken by considering 
higher order moments of the velocity distribution function. 
Though not shown, we have also explored the effects of vary- 
ing the specific dark halo profile shape by varying a in equa- 
tion 2. As a increases, the allowed core size increases at fixed 
Vmax- For example, at Vmax = 30 kms"' the best-fit value of 
'"core varies from 50 pc to 250 pc as a varies from 1.0 to 2.0. 
We discuss the origin of this variation in the next section. 

Note that while some large values of rcore (^ kpc) are al- 
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lowed by the velocity dispersion data, they require very large 
values of Vmax 100 kms"'). One may be concerned that 
the quoted Vmax values represent an unphysical extrapola- 
tion to large rmax radii based on constraints derived by stel- 
lar velocity dispersions at ~ 2 kpc. However, direct obser- 
vational evidence suggests that this is a reasonable extrapo- 
lation. Not only is our adopted profile shape motivated by 
rotation curve fits (Burkert 1995) but large dark matter ha- 
los around isolated galaxies are known to decline no faster 
than at large r (Prada et al. 2003). In this case, the only 
effect that could act to truncate this fall-off is tidal strip- 
ing. If Fornax is at a distance D from the Milky Way, then 
the tidal radius may be estimated via the Jacobi approxi- 
mation /?, ~ Z) (MFornax« R,)/i3MMwi< D))) For a flat. 
Milky Way rotation curve with Vmax = 220 kms~' we obtain 
R, ~ 80 kpc (Vmax/(200 kms-l))^^^ Compare this to our cen- 
tral Vmax — 200 kms~' solution, which gives rmax = 3.20 = 
32rcore — 48 kpc and we see that the peak rotation curve 
should be little affected by tidal truncation. We conclude 
that while a large Fornax halo may be implausible on physical 
grounds (see below) the derived constraint itself is meaning- 
ful. 

It is also possible that the large Vmax end of the contour 
is affected by the presence of interloping stars, and that the 
contours actually close at a smaller value of Vmax (and rcore)- 
In particular, removal of interlopers is most important in the 
outer regions of the galaxy. To examine the effect of inter- 
lopers, we have performed an independent fit of the data pre- 
sented in Walker et al. (2005). We use the same number of 
radial bins, and in each bin we assume that the distribution 
of stars in velocity follows a Gaussian plus a constant, where 
the constant accounts for interlopers. The main difference we 
find is that there is a lower dispersion in the last bin at ^ 2 
kpc than in Walker et al. (2005), with a larger error in this 
bin. These results do not change the shape of the constraints 
at low Vmax, but disfavor large values: Vmax ^ 65 kms~^ and 
''core ?J 400 pc are excluded. 

While uncertainties associated with interlopers stars may 
affect the closure of the contour in Figure 2 at large values 
of Vmax, the lower part of the contour is more robust to these 
uncertainties. Indeed, it is this region of the diagram, and the 
Vmax-'"core degeneracy itself which allows us to use the implied 
phase-space densities in these galaxies to limit core sizes. Be- 
fore going on to discuss these limits, we discuss the origin of 
this degeneracy using analytic arguments and extend these ar- 
guments to explore the dependency of our constraints on the 
assumed underlying dark profile shape. Readers interested in 
our resulting constraints on the core size of Fornax may skip 
ahead to section 3. 



2.1. Toy Model to Understand the Solutions from Stellar 

Kinematics 

Here we show that the scaling of the contours in the Vmax- 
''core can be understood in the context of a toy model. For 
simplicity we will consider /? = 0. From equation 1 we see 
that if the stellar mass is sub-dominant, the stellar velocity 
dispersion is governed by the circular velocity curve of the 
dark matter Vdr) within the stellar radius, r^ ^ rc ^ 0.5 kpc. 
Here r^, corresponds to a weighted average of the stellar mate- 
rial over an integral / drp*(r)Vj?(r)/r. Thus r,^ will typically 
be smaller than the half-mass radius of the stars. 

Consider now a toy-model dark matter halo with a rotation 



curve given by 



Vc(r) = 



flVm 



flVm 

v„„ 



r < r., 



y 'core y 



1/2 



<r< rflat 

r > Tflat = Tcore/fl^ 



(5) 



This corresponds to a density profile with a constant density 
at r < fcore- The profile transitions to an p oc r"' regime be- 
yond the core and finaUy reaches an isothermal distribution, 
p oc r"^, at r > rflat. We assume that the transitions occur in a 
small region around rcore and r^-^i and that the density profile 
is continuous. The variable a < I parameterizes the sharp- 
ness of the transition from a core to flat rotation, and can be 
viewed as analogous to the parameter a in equation 2. For 
a = 1, the transition is sudden, rcore = ''flat, while for a ^ 1 
the rotation curve rises very gradually and reaches Vmax at 
r = rant = '"core/a^- Note that the value of the rotation veloc- 
ity atr = rcore is Vcok = «Vmax, and thus decreases at fixed Vmax 
as the transition is made more gradual. Also note that the den- 
sity in the core of the halo po is proportional to a Vmax /''core- 

Now consider a stellar distribution with radius r^, and 
with constant density, p*(r) =constant, embedded within the 
gravitationally-dominant dark halo. If the halo has a very 
sharp transition from the core region to flat rotation, a = 1 
and rcore ' ''flat, then equation 1 implies that the central steUar 
radial velocity dispersion, = cr,t(r = 0), is 



I ^max 



(w) 



r* < rc, 



1 11 

Vmax [\ +ln(r^/rco,«)] > rcore = '•flat- 



(6) 



This suggests that if we fix the central velocity dispersion, the 
maximum circular velocity of the halo is degenerate with the 
core radius, and must increase with core radius as Vmax oc rcore 
in the large-core regime. Another way to interpret the result 
is that in this regime, the core density of the dark matter halo 
is fixed, but not the core size. Conversely, in the small core 
regime we expect Vmax to vary slowly with increasing core 
size ( Vmax oc ^1 -1-2 ln(r-,/ rcore) for this example). Interest- 
ingly, these qualitative features are demonstrated by our best- 
fit contours in the right panel of Figure 2. 

For a more gradual transition from core to flat rotation, we 
allow a < 1 with rcore < ''flat- In this case, for r^ < rflat, we 
have 



^^Vinax ( r ) ^ ^ core 



aVm 



'* core 2 



1/2 



(7) 



> ra 



This means that for a fixed core size, a weaker transition 
(smaller a) requires a larger Vmax to reproduce a given a^: 
Vmax OC rcore/ a or fixed halo core density in the large-core 
regime, and Vmax oc ^JT^/a in the small-core regime (rcore ^ 
r*). This result demonstrates qualitatively why the implied 
core sizes increase with a sharper transition in the assumed 
dark matter profile (a : 1 — > 2 in equation 2). The stellar ve- 
locity dispersion is roughly governed by Vdr ~ r*). As the 
transition becomes sharper, the ratio Vdr = r^)/Vm!a increases, 
and an observed steUar velocity dispersion can be produced 
with a smaller Vmax- 

The scalings seen in the left-hand panel of Figure 2 for g = 
Po/'^dm Vmax cau be similarly understood. In the large- 
core regime, where r^, < rcore, the central density po will be 
approximately constant because fixing roughly fixes Vdr*)- 
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In this case, we expect Q oc a~^^ oc (aV'max)"'' « r~^^^, which 
reproduces the contour scaHng in Figure 2. 

The dependence of Q on the transition parameter a is 
weaker than the scaling of Vmax: cx; aymax''core[l ~ 

21n(fl)]-'/^. This result is consistent with our numerical find- 
ings: the sharper the transition, the larger is the core size re- 
quired for a given Q. In the small core regime, we have seen 
thataVmax 0^ \/^core and hence Q oc ^^ore implyiiiS the fcort 
vs Q contour should steepen slightly for small values of rcore- 

3. RESULTS 

3.1. Phase-Space Constraints 

The derived relationship between allowed Vmax and rcoie val- 
ues becomes interesting when expressed in terms of the im- 
plied central phase-space density for each solution. Consider 
the case of WDM with a distribution function given by 

where g is the number of spin degrees of freedom of the par- 
ticle. For /3 = 1 the distribution is thermal, and in this case we 
define the mass of the WDM particle as m^. An example of a 
(non-thermal) /3 ^ 1 distribution is WDM from an oscillation- 
produced sterile neutrino (e.g. Dodelson & Widrow 1994). If 
/3 is independent of momentum, we may calculate the max- 
imum phase-space density by fixing the present-day density 
of dark matter and calculating the velocity dispersion from 
equation 8: 

emax = 5 X 10-V (I) (^-^)'M0pc-3(kms-')-l (9) 

As discussed in the introduction, the (course-grained) phase- 
space density in any WDM galaxy halo can never exceed 

Gmax- 

Measurements of the Lya power spectrum place a limit 
on the free streaming scale (or equivalently, the present-day 
velocity) of WDM particles. We can convert this directly 
to a lower limit on Q. For a fixed free-streaming scale, 
the mass of a thermal and non-thermal particle is related by 
fi^x/Tx = nts/Ts, where mj is the mass of the non-thermal par- 
ticle and Ts is its corresponding temperature in equation 8. 
Fixing the present-day density of dark matter in this case im- 
plies f3 = (ms/nijc)* as long as (3 is independent of momen- 
tum. Note that since P and Q both scale as the fourth power 
of the WDM particle mass, any Lya forest constraint can be 
mapped uniquely to a Qmax constraint, independent of the par- 
ticle mass. 

The most conservative quoted constraint from the Lya for- 
est power spectrum expressed as a limit on the thermal WDM 
particle mass is nix ^ 0-5 keV (Viel et al. 2005; Abazajian 
2006b). A distinct analysis using different simulations has 
derived an even stronger limit > 2.5 keV (Seljak et al. 
2006). The > 0.5 keV limit implies an upper limit on 
the phase-space density in the core of any galaxy of Q > 
3 X lO~5M0pc-^(kms~')~^ That is, if WDM accounts for the 
mass of the dark halo of Fornax, then Figure 2 implies that the 
core must be quite small, rcore ^ 85 pc, and Vmax ^ 35 kms"' 
(the lower dotted line). The more stringent 5c; 2.5 keV 
constraint implies Q > 0.02 M0pc"^(kms"')~^ and demands 
''core ^ 10 pc even with allowances for extremely sharp tran- 
sitions to the profile core. 

Of course, the Lya forest constraint only applies to the spe- 
cific case of a WDM particle with a phase- space density given 



by equation 9. More generally, there are models with a fi- 
nite primordial phase-space density that do not have the same 
mapping between Q and Pik) as does WDM. For example, the 
well-motivated SuperWIMP scenario suggests that dark mat- 
ter arises from out-of-equUibrium decays and this gives rise to 
a non-thermal spectrum (Cembranos et al. 2005; Kaphnghat 
2005). Also, non-thermal resonantly-produced sterile neu- 
trino dark matter would be "cooler" than a thermal or non- 
resonant sterile WDM candidate (Shi & Fuller 1999; Abaza- 
jian et al. 2001). Yet another example is fuzzy dark matter 
(Hu et al. 2000). 

In cases such as these the relation between the phase-space 
density and the power spectrum is different from that im- 
plied by equation 9, a more direct observational constraint 
on Q may be derived from considering the rotation curves 
of low-mass spiral galaxies (Dalcanton & Hogan 2001; Ko- 
rmendy & Freeman 2004). Among the highest-resolution 
rotation curves for low-mass spirals are those presented by 
Simon et al. (2005), who used two-dimensional Ha veloc- 
ity fields to extract rotation curves for five galaxies. When 
the dark matter halo components of these galaxies are fit to 
cored density distributions the central densities range from 
po ~ 0.1-0.5 Mqpc"^ and the rotation curves flatten at 
Vmax ^74-114 kms~'. If we adopt adm — 0.55 Vmax (as is 
appropriate for the cored profiles we have used here) then the 
implied central phase-space limits are 2 > (1 .7, 4.4, 1 .7, 1 .6, 
and 0.5) x lO"*^ M0pc-3(kms-i)-3 for NGC2976, NGC4605, 
NGC5949, NGC5963, and NGC6689, respectively. If we 
(conservatively) adopt Q > 10"^ M0pc~^(kms~')~^ as a (di- 
rect) constraint on Q, this provides a fairly clean upper-limit 
on the Fornax core radius, rcore ^ 300 pc, for a broad class of 
dark matter candidates (upper dotted fine in Figure 2). 

3.2. Phase-Space Independent Constraints 

While the most popular alternatives to CDM produce cores 
because of their phase-space constraints, cores in dark matter 
halos need not arise as a result of large velocity dispersions. 
For example, self-interaction or self-annihilation would pro- 
duce cores even if the dark matter is cold (Spergel & Stein- 
hardt 2000; Medvedev 2000; Kaplinghat et al. 2000; Sanchez- 
Salcedo 2003). Astrophysical or dynamical effects could also 
produce cored profiles. For our most general constraint we 
can use the rcore-Vmax relation demanded by the velocity dis- 
persion data to place a reasonable upper limit on the core size 
in Fornax by imposing a plausible upper Umit on Vmax for its 
halo. 

Models of Fornax with large Vmax have large masses. As 
a result, they will experience large dynamical friction and 
spiral to the galactic center very fast. Therefore only in the 
rare circumstance of a very recent accretion event could a 
massive subhalo at a distance of ~ 140 kpc from the Milky 
Way be possible. Zentner & Bullock (2003) found that only 
5% of Milky Way-sized ACDM halos contain a subhalo with 
Vmax > 100 kms~' and that the fraction falls sharply beyond 
this point. Since the presence of dark matter cores would 
only decrease the subhalo population, we may adopt Vmax ^ 
100 kms"' as a reasonable upper limit on the Fornax halo. 
This gives rcore 700 pc as a conservative general upper limit 
on the core size of Fornax. 

In Ught of this argument, a core large enough to explain the 
extended distribution of globular clusters in Fornax, rcore — 
1.5 kpc, seems highly implausible. Such a case would de- 
mand a very large Fornax dark halo: Vmax 200 kms"' with 
Mi< rmax) — 4.4 x 1O"M0. The dynamical friction timescale 
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for such an object at a distance of ^ 140 kpc is less than 
^ O.SGyr. Consider also the fact that the stellar mass of For- 
nax is ~ 3 x IO^Mq (Mateo 1998). This means that a 
1.5 kpc core in Fornax would require both that the Milky 
Way is experiencing a (very rare) recent equal-mass merger 
and that the merging object is a system with a stellar baryon 
fraction that is ~ 10,000 times smaller than the universal 
baryon fraction. 

4. CONCLUSIONS AND DISCUSSION 

In this paper, we have shown that a standard WDM particle 
cannot produce a core in Fornax larger than ~ 85 pc without 
saturating the impUed Lya forest bound on the phase-space 
density. For a more general class of dark matter, the rela- 
tion between the phase-space density and the power spectrum 
may be different, and in these cases the limit on Q comes from 
low-mass galaxy rotation curves, giving rcore 1^ 300 pc. Irre- 
spective of the origin of the core, dynamical arguments alone 
place an upper limit on the Fornax halo mass and demand 
W < 700 pc. 

One important caveat to the above results is that there is no 
concrete prediction for the precise shape or normalization of 
dark matter halo profiles in cases where the dark matter is not 
cold. Even in well-studied cases like WDM, a phase-space 
core has never been realized in an N-body simulation and the 
nature of the expected transition between a rising p{r) profile 
and the flat core region is unknown. The density profile we 
have adopted in equation 2 asymptotes to the CDM predic- 
tion at large radius, p{r) oc and allows some freedom in 
the nature of the transition from the inner core region to the 
outer fall-off in density. The sharpness of this transition is 
captured in the shape parameter a: the sharper the transition 
(larger a), the closer the core radius is to the radius of maxi- 
mum rotation. Given our general ignorance as to the "correct" 
shape, we have explored the degree to which our constraints 
on the core size change when we adopt a = 1 and 2. The 
Q constraint for a broad class of dark matter models derived 
from from low-mass spiral galaxies yields rcore ^ 500 pc and 
200 pc for a = 2 and 1, respectively. For WDM, the Lya for- 
est power spectrum implies rcore 150 pc and 40 pc, for a = 2 
and 1, respectively. 

What do our results mean for the globular clusters in For- 
nax? A core large enough to explain the extended distribution 
of globular clusters, rcore ^1-5 kpc, would imply a dynami- 
cally implausible dark halo, Vmax ^ 200 kms~' (see Figure 2). 
In addition, such a core is vastly inconsistent with any viable 
WDM model. 

Given these constraints, an alternative explanation for the 
globular cluster system may be considered. For example, the 
profile could be cuspy, but this would require the globular 
clusters to have formed at large radii beyond the asymptotic 
central slope (J. Read, private communication). This seems 
to require a fine-tuned timing argument. The timing argument 
could be sightly alleviated by arguing that only the irmermost 
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globular cluster defines the core (rcore ^ 240 pc), however 
even this case is inconsistent with standard WDM. 

We may also consider a solution based on tidal heating. For 
example, dynamical friction effects on the orbital decay of 
the globular clusters can be counteracted by dynamical heat- 
ing via tidal shocks or by relaxation during major mergers 
(Oh et al. 2000). The tidal heating explanation may have dif- 
ficulty because the proper motion measurements indicate that 
the Fornax is currently near its orbital pericenter and has apoc- 
enter of ~ 250 kpc (Dinescu et al. 2004) and thus does not 
come sufficiently close to the center of the Milky Way to ex- 
perience strong tidal force. Note, however, that the tidal heat- 
ing could have been provided by another object in the past 
and not necessarily by the Milky Way progenitor (see, e.g., 
Kravtsov et al. 2004). Heating by a major merger with another 
dwarf galaxy is plausible, as there is evidence for such recent 
major merger in Fornax from observations of extra-tidal stars 
(Coleman et al. 2005; Olszewski et al. 2006). Although such 
mergers are not likely for satellites orbiting in a potential of a 
much larger host halo, given large pericentric and apocentric 
distances of the Fornax, it is likely that it has been accreted by 
the Milky Way only recently. Major mergers are much more 
hkely prior to accretion by a large host. Other possible, albeit 
more speculative, explanations, such as core formation due to 
inspiralling of multiple black holes, are discussed in Sanchez- 
Salcedo et al. (2006). 

As emphasized by Goerdt et al. (2006) and Sanchez- 
Salcedo et al. (2006), the Fornax globular cluster system 
presents an interesting case study with possible fundamen- 
tal implications for the nature of dark matter. Here we argue 
that interpretation in terms of large dark matter core is prob- 
lematic, as such a large core is not consistent with popular 
dark matter alternatives to CDM. Large cores also appear to 
be inconsistent with high-resolution measurements of rotation 
curves of the LSB galaxies for an even broader class of dark 
matter scenarios. Dynamical heating of the globular clusters 
by tides or relaxation in mergers could provide an alternative 
explanation and needs to be explored. 
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